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Abstract 



In the previous paper we studied Dirac fermions in a non-Abelian random vector 



potential by using lattice supersymmetry. By the lattice regularization, the system 



of disordered Dirac fermions is defined without any ambiguities. We showed there 
> 

that at strong-disorder limit correlation function of the fermion local density of states 



decays algebraically at the band center. In this paper, we shall reexamine the multi- 
flavor or multi-species case rather in detail and argue that the correlator at the band 
center decays exponentially for the case of a large number of flavors. This means 
that a delocalization-localization phase transition occurs as the number of flavors is 
increased. This discussion is supported by the recent numerical studies on multi- 
flavor QCD at the strong-coupling limit, which shows that the phase structure of 



QCD drastically changes depending on the number of flavors. The above behaviour 



/^ • of the correlator of the random Dirac fermions is closely related with how the chiral 

symmetry is realized in QCD. 
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1 Introduction 

In the last several years, Dirac fermions in a random vector potential and/or with ran- 
dom mass have been studied by various methods. This system is closed related with 
the phase transition between plateaus in the quantum Hall states and quasiexcitations 
in the d-wave supercondctor|], §, |], H; H H; 0; H- However it is still controversial if 
there exist extended states at the band center. In the previous paper 0, we studied 
the case of non-Abelian U(N C ) random vector potential for large N c . There both 
weak and strong disorder cases are investigated separately by using lattice super- 
symmetry(LSUSY) which was originally invented for the study on the SUSY gauge 
theory. Especially it was shown that correlation function of the fermion local den- 
sity of states(DoS) decays algebraically at the band center. This result indicates the 
existence of extended states there. 

The above random Dirac fermion system on the lattice is closely related with 
the lattice QCD for which there exist many useful and definite studies by numerical 
simulations. Among them, one of the most interesting studies shows that the phase 
structure of lattice QCD changes drastically depending on the number of flavors of 
quarks HlO|. For small number of flavors, quarks are confined and the chiral symmetry 
is spontaneously broken, whereas as the number of flavors is increased, quarks are 
deconfined and the chiral symmetry is restored. Then mass of "pions" is nonvanishing 
even for the vanishing bare quark mass for many flavors. 

The above result of multi-flavor lattice QCD holds even at the strong-coupling 
limit which directly corresponds to the strong-disorder limit in the random Dirac 
fermions. In the previous study on the Dirac fermions in a random U(N C ) vector 
potential 0, we implicitly assumed a small number of fermion flavors. In this paper 
we shall reexamine the system of multi-flavor or multi-species Dirac fermions. The 
flavor degrees of freedom is sometimes introduced in the replica tricks and one takes 
the limit iVy(=the number of flavors) — > after calculation. However in the present 
study, Nf is kept finite and moreover we even consider the case like Nf/N c > 1. In 



the above lattice QCD, N c = 3 and the critical value of Nf is 7. 

Plan of the present paper is as follows. In Sec. 2, we shall study the strong- 
coupling limit of lattice QCD by deriving an effective action of U(N C ) invariant 
excitations, which is a kind of antiferromagnetic spin model. From this effective 
action, we can study the vacuum structure, especially how the chiral symmetry is 
realized. For a small number of flavors, the chiral symmetry is spontaneously bro- 
ken and (quasi) massless pions appear as a Nambu-Goldstone boson as we showed 
previously |11||. For a large number of flavors, on the other hand, we show that the 
chiral symmetry is restored and "pions" become massive even at the chiral-invariant 
limit. This result is of course in greement with the numerical studies^]. In Sec. 3, 



we review the lattice model of the Dirac fermions in a U(N C ) random vector poten- 
tial. Especially we explain the LSUSY which is introduced in order to take ensemble 
average over the random vector potential. Then we shall derive an effective action 
by using similar techniques used for QCD. From the discussion on QCD in Sec. 2, we 
show that there is a phase transition as the number of flavors is increased. This phase 
transition can be observed from the correlator of the fermion local DoS. Section 4 is 
devoted to conclusion. 

2 QCD with many flavors and chiral symmetry 
restoration 

Before going into detailes of the Dirac fermions in a random vector potential, we 
shall investigate the strong-coupling limit of the lattice QCD. We follow the dis- 
ussion in Ref.jnj which is also applicable for the random fermion system rather 
straightforwardly . 

We consider d- dimensional hypercubic lattice. Action of the QCD at the strong- 



coupling limit is given by 



Sqcd = - Y, [i){x)-i^U^{x)i){x + /i) - ip(x + ^^Ul{x)i)(x) 



(2.1) 



where x = (x®, ..,Xd-i) denotes lattice site, /i = (0, ...,d — 1) is the direction index, 
U^{x) is U(N C ) field on the link (x, x + jj) U^{x) = (L^(x)J £ U(N C ) and we set the 
lattice spacing at — 1. The Dirac fermion ip carries not only the spinor and color 
indices s = 1, ..., 2^ d ' 2 \ a — 1, ..., N c but also flavor index 1 = 1,..., Nf, i.e., ip = "0"' 7 . 
In Eq. fl2.ip , 7 M is the 7- matrices which work on the spinor index of ip. We add the 
following quark mass term in the action 

S M = -M B J2iJ>( x M x )> ( 2 - 2 ) 

where Mb is the quark bare mass. The gauge field action is vanishing because we 
are considering the strong-coupling limit. In the numerical studies [10|, the Wilson 
fermion is used for defining definite number of flavor. Here we consider the naive 
lattice fermion. However it is quite plausible that a similar phase transition occurs 
as "the number of flavor" Nf is varied. 

It is well-known that Sqcd i n (|2.1|) can be rewritten more compact form by the 
following change of variables, 

ip(x) = T(x)x(x), 4>(x) = x(x)T f (x), (2.3) 

with T(x) = ( r yo) x ° ■ ■ ■ {'jd-i) Xd ~ 1 and using the identities like (7 M ) 2 = 1 and (7o) n 7i = 
( — ) n 7i(7o) n (n is an integer), etc, 

Sqcd = -= J2 [x( x )v»(x)U^(x)x(x + fi) - x(x + /i)r7 M (x)£/J(x)x(x) 

= J2x{x)D X (x), (2.4) 

where rj (x) = 1, rji(x) = (— ) x °, r] 2 (x) = (— ) Xo+Xl , ■ ■ ■, and 

Dx(x) = - Y, [v^(x)U^{x)x{x + fi)- Vn( x - v) u t(x - fi)x(x - //)] . (2.5) 



Similarly 

S M = -M B J2x(x)x(x). (2.6) 

In Eqs. (|2.4j) and ( |2.6| ), the spinor indices of x an d X are diagonal and they play the 
same role as the flavor index. 

Partition function of the system is given by 

Zqcd = J[DUDxD X ]e s ^ +s ^, (2.7) 



where [DU] = TlimkdU^x) is the Haar measure of U(N C ). In Eq.( |2.7| ), the integral 
over the gauge field U^(x) can be performed by the 1/N C expansion because it reduces 
to the one- link integral in the strong-gauge-coupling limit [[11]. There are two "phases" 



in the one- link integral 12]. The present case corresponds to the "strong-coupling 
regime" . Detailes of the calculation can be seen in Ref . || [TT[] and here we simply give 
the final expression. To this end, we define the following color singlet quark bilinear 
operators, 



1 



m /)W 



TrEx a '»W(*) 5 (2-8) 



Nc a 



where a and (3 are flavor-spinor indices which take Nf s = Nf x 2^^ values. It is 
obvious that after integral over U^(x), the partition function Zqcd is expressed in 
terms of m,p(x), 

Zqcd(J) = J[D X Dx]e s ^ +SM+J - m , (2.9) 

where J's are source fields and J ■ m = Tr( .Ira) = J2 Ja' m a- m ^ ne leading-order of 

1/iVeQ 

i-S m (A) = £TV[y(A M (x))], (2.10) 

with 






Xfi^x) = m(x)m(x + jj), 

g(X) = 1 - (1 - A) 1 / 2 + In [~(1 + (1 - A) 1 / 2 )] , (2.11) 



lr rhere is a systematic expansion in powers of 1/N C . However, essential feature of the effective 
action and the final results are not influenced by the higher-order terms. 



and Tr is the trace over the flavor-spinor indices.^ 

The path- integral over x an d X in Eq. (|2.9|) can be performed by introducing ele- 
mentary meson fields Mg(x). In terms of Aip(x), the partition function is expressed 

as 

Z Q cd(J) = J[DM]e s *ffM +J - M , 
^S eff (M) = Y,^[9Mx))]-'E[TrlnM(x)-M B T r M(x)l (2.12) 

where A M (x) = A4(x)A4(x + /x) and integral over M. is defined by polar decomposing 
A4 as M. = RV with positive-definite Hermitian matrix R and U(Nf s ) matrix V, 
and JdM = J dV (=the U(N fs ) Haar measure) [0] . 

From the effective action ([2.121) , we can study the vacuum structure of the system [ITT 



"Expectation value" of M. (x) is obtained by the effective potential which is obtained 
by substituting M.p(x) = vSq in S e ff, 

Ve f f(v 2 ) = -dg{v 2 ) + ~ hit; 2 - M B v. (2.13) 



N c N fs ^ ' ^ ' 2 

From (|2~T3| ), 

v = f^ + 0(M B ). (2.14) 

For vanishing bare quark mass M B = 0, S e ff(M) is a function of A M (x), and then 
the vacuum expectation value of the fields M. (x) is given as 

, vV at even site 
(M(x)) -- { (2.15) 

vVq at odd site, 

where Vo G U(Nj s ). From the above discussion, low-energy excitations are obtained 
by introducing "pions" kk{x) as 

, vV(x) = ve i ^-<K' KK( - x ^ TK at even site 
M(x) = { K J (2.16) 

vV\x) = ve-^K^Kix)^ at odd gite5 



2 We often omit the flavor-spinor indices on m(x) etc. 



where T K, s are generators of U(Nf s ). 

For N c /Nf > 1, we can expect that fluctuation of ttk{x) is not large as S e ff(JA) oc 
N c and topologically nontrivial configurations do not contribute to the partition func- 
tion. Therefore we expand V(x) as 

V(x) = 1 + z^k(x)T k + ---, (2.17) 

K 

and then 

^S eff (M) ~ V^E(V^W) 2 - M B vJ2^ K (x). (2.18) 

Interaction between the pions is suppressed by 1/N C . In this case, the chiral symmetry 
is spontaneously broken and 7Tk-(x)'s are quasi-Nambu-Goldstone bosons [pi] . 



What happends for N c /Nf s <C 1. Numerical studies indicate that the "pions" 
become massive even in the limit Mb —* 0. This means that chiral symmetry is 
restored in this case. Actually this is not so difficult to see this behavior in many- 
flavor case from S e ff. The effective action ( |2.12| ) can be regarded as a system of 



U(Nf s ) antiferromagnets.Q One can expect that as degrees of freedom of the U(Nf s ) 
spin, i.e., V^ in the present system, becomes large, the long-range order is destroyed 
by fluctuations. For example in spatial 2 dimensions, ordered phase exists in the Ising 
model whereas there is no ordered phase in the 0(3) spin model. Therefore we can 
naturally expect that the chiral symmetry is restored for Nf s /N c ^> 1. Without the 
spontaneours breaking of the chiral symmetry, the expansion (|2.17| ), which assumes 
small fluctuation of ttk(x), is not justified. 

Let us see the above chiral phase transition more closely. To this end, we shall 
introduce a more tractable model which is expected to belong to the same universality 
class of S e ff. Let us notice that as M B — > 0, S e ff(M) is a sum of the one-link term 
like 



Hink = N c Tr g{\^x)) - — In X^x) . (2.19) 



1 

2d 

3 More direct relationship between the lattice QCD at strong coupling and antiferromagnets can 

be seen in the Hamiltonian formalism]! 31. 
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Then we consider the following model which consists of a similar one-link term, 



S v = — J2 Tr\V*(x)V(x + n) + H.c. 

90 x,n 



(2.20) 



where the parameter — ~ N c v 2 . Actually (a continuum version of) the action 
Sy is derived directly from the QCD action ( j2.1| ) by using Zirnbauer's color-flavor 
transformation |TJ| (see also Ref.[15|]), and therefore it is quite plausible that S e ff and 
Sv belong to the same universality class. 

The system Sv ( |2.20|) can be easily studied by the mean- field theory(MFT) . To 
this end, let us assume the expectation value of the field V(x) G U(Nf s ) as 



(V^(x))=w6l 



(2.21) 



where we have used the U(Nf s ) symmetry and we assume that w is real.[] We employ 
the following one-site action as MF action |I6[, 



Smf — — Tr 
9o 



V + V ] ] N fs w 



9o 



(2.22) 



and the partition function of the MFT is given by 



J MF 



-V M f{w) 



dVe 



Smf 



(2.23) 



From Eqs.( |2.22| ) and ( |2.23| ), it is obvious that the stationary condition of Vmf{w) 
gives the "gap" equation which determines the expectation value w. 

The integral over V in ( |2.23| ) with the U(Nf s ) Haar measure can be performed by 
the same techniques used before for integrating over the gauge field U^x), 



VWM = _*,.,(«£)>) + -V..> 



Especially for small w, 



V M f(w) 



*Nfs go' 



N f 



w 



9o 



— — H —w+Olw). 

\9oJ 9o 



(2.24) 



(2.25) 



4 Please do not confuse w with v. In the terminology of the spin model, v is the magnitude of 
spin, whereas w is the magnetization. Therefore nonvanishing w means the long-range order. 



From Eqs.( |2.24] ) and ( |2.25| ), it is obvious that for Nf s > — = N c v 2 , w = is the stable 



minimum, whereas for Nf s < — = N c v 2 , w = becomes unstable and w has a finite 
value at the minimum of Vmf- Then this simple consideration using Sy and the MFT 
supports our expectation, i.e., for Nf s /N c <C 1 the vacuum is chiral symmetric and 
there are no Nambu-Goldstone bosons, whereas for N c /Nf s <C 1 the chiral symmetry 
is spontaneously broken and gapless excitations exist. 

In the chiral-symmetric case, physical quantities are calculated by using the U(Nf s ) 
Haar measure for V(x). We are interested in the meson-meson correlation functions 
like {Aip(x)A42,(0)) . The integral over V(x) can be performed by using the character- 
expansion for V(x)V\x + jj) G L/"(iYf s )P7|) and each link term of the action can be 
written as; 

exp(S /mfc ) =Y.F p (N c ,v) Xp (V\x)V(x + ri), (2.26) 

p 

where p refers to irreducible representations of U(Nf s ), xp is the character and 
F P (N C , v) is a function of N c and v. 

Integral over V G U(Nf s ) can be performed by the following formula [|H]j. 

rfrVJ = -^5 a J}. (2.27) 

From Eqs. (|2.26|) and ( 2.27 ), the meson correlation function is calculated for the sim- 



/ 



plest case, i.e., one-dimensional system with free boundary condition. In this case the 
V^V correlator is obtained as follows; 

= e~ m ^ 8Z 6%, (2.28) 

where / and refer to the fundamental and trivial representations of U(Nf s ), re- 
spectively, and m n = ln(Nf s Fo/Ff). As the dimension of the group U(Nf s ) becomes 
larger, the mass m n is getting larger as we expected. 

The above considerations naturally explain the numerical calculations of the multi- 
flavor QCD[|lO|. In the following section, we shall apply the above observation to the 
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Dirac ferimons in a random vector potentail. We shall begin with reviewing the 
previous studies using LSUSY|J. 

3 LSUSY and localization 

Let us start with reviewing lattice formulation of random Dirac fermions. Our model 
is motivated by the network model for the quantum Hall phase transition in two- 



dimensional electron systems in a strong magnetic field [19]. There electrons move 
along equipotential line of the random impurity potential and acquire Aharonov- 
Bohm phase from the strong magnetic field. Then we shall consider d = 2 case in 
this section. Fermionic (electronic) part of the action is given as in the lattice QCD 
in the previous section; 

S D = J2x(x)D X (x), (3.1) 

where xi x ) an d xi x ) are two-component Grassmann fields with the color U(N C ) and 
flavor U(Nf s ) indices as before. The bare mass Mb in Sec. 2 corresponds to the energy 
measured from the band center in the present electron system and we set Mb = 
as we consider states at the band center. Corresponding to the Grassmann fields, we 
introduce complex scalar fields 4>(x) and <f>*(x) with the color and flavor degrees of 
freedom as follows [jlS|[ , 

S^Y.D^^D^x). (3.2) 

It can be shown that the system So + S^ possesses LSUSY, 

6(f) = e + x- + e-X+, 

b<t> ] = X-£+ + X+e-, 

$X± = -D<pe±, 

5 X ± = -Dtfe ± , (3.3) 

where x± is the chiral decomposition, i.e., 7sX± = ±X± an d e± are anticommuting 
spinor variables with chirality 7 5 e± = ±e±. From the LSUSY ( |3.3|) , the partition 
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function is just unity, i.e., the determinant from the x integral and that from 
integral cancell with each other |§. Parameter g which controls randomness of the 
vector potential can be introduced as follows, 

riV j. i 

P[U] = []exp p tr(U,(x) + 0j(x))]. (3.4) 

x 9 

Ensemble average with respect to the random vector potential for calculating corre- 
lators etc can be performed readily by the LSUSY||; 

(O) = [[DUDxD X D^D(j)]P[U]e- SD - s *0 

lDUD ^ ]m mm] ° 

= JiDUD^mul^^O, (3.5) 

where O is physical observable of the Dirac fermion and |%; U) is the lowest-energy 
state of x( x ) i n the background of U^(x). From Eq. (|3.5|) , the utility of the LSUSY is 
obvious. In the rest of the discussion, we consider the strong-disorder limit g — > oo 
and therefore there is no distribution weight for the gauge field U^x) as in the strong- 
coupling limit of QCD. 



The bosonic action Sa can be rewritten into more compact form 1 18 



S^ =>S UV = J2 {^D<p + (p*Du] , (3.6) 



where u and tp are complex boson fields. Then the one-link integral of U^x) in 
the previous section can be applied straightforwardly. Fermion-fermion(FF) sector is 
given by S m (X) in Eq.( p.lOj ). Boson-boson(BB) sector is also written in terms of the 
following composite fields like jj- Yl a , { P a,a { x )'ra8{ x )i ^- Yja u ° a,a { x ) UJ a,p{ x )-> e ^ c --> but 
the sign of the action is just opposite to the FF sector. There are additional terms 
which include fermionic composites like x a ( x )f a ( x ) but they contributions to the FF 
and BB sectors are negligible in the leading order of 1/N C . 

From the above arguments and investigation on QCD in Sec. 2, it is straightforward 
to discuss physical properties of the random Dirac fields at the band center which 
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corresponds to the vanishing bare mass Mb — > 0. We are interested in the correlation 
function of the local fermion DoS m^(x) = J^a^ai^^ix), which is nothing but 
the "pion" field in QCD. Therefore in the case of a small number of flavors, one 
expects that the operator m^(x) condenses and has a nonvanishing expectation value. 
However as we are studying the 2-dimensional system, the genuine long-range order 
does not exist and the correlator of m^(x) decays algebraically (quasi-long-range 
order) instead |§|] 

(^(x)^(0)) ~ \x\~ 1/Nc . (3.7) 

The power-decaying behavior in (|3.7|) comes from the fluctuation of the "pions" n %:{%)- 



It should be remarked here that in the present lattice formulation topologically non- 
trivial configurations of the random vector potential are all included. This formulation 
is motivated by the network model by Chalker and Coddingtonjnj and it is in con- 
trast to other formulations like that using the conformal field theory. Therefore it is 
not so surprising even if the decay power of the correlator ( |3.7| ) is different from that 
obtained in other methods^]. The power-decay behaviour ( |3.7| ) means the existence 
of the extended states anyway. 

Let us turn to the many flavors case. The discussion in Sec. 2 can be applied 
straightforwardly to the random system in this section. As the system is 2-dimensinal, 
the disordered phase is enhanced in the U(Nf s ) antiferromagnets. The correlator of 
the local fermion DoS behaves as 

(^(x)'ipip(O)) ~ e~ MM , (3.8) 

though we cannot estimate the "inverse localization length" Mi analytically. The 
above is a new and main result in this paper which indicates that extended states at 
the band center disappear as the number of flavors is increased. 



Recently certain related model is studied in Ref. [20|, and a similar phase transition 



is predicted, though the model is defined in the continuum space instead of on the 



D Here we do not show explicitly the flavor-spinor indices as they are obvious as in (2.28) 
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lattice. In the present lattice model, which is motivated by the network model for the 
quantum Hall transition, topological nontrivial configurations of the gauge field are 
included. Therefore it is quite nontrivial that a similar transition or crossover occurs 
in the both cases. 

4 Conclusion 

In this paper we studied Dirac fermions in the random U(N C ) vector potential by 
using lattice formulation. Especially we showed that physical properties at the band 
center drastically change as the number of flavors is increased. For a small number 
of flavors, the correlator of the local fermion DoS exhibits algebraic decay, whereas 
for many flavors, it decays exponentially. This result suggests that extended states 
at the band center disappear as Nf — darge. Numerical studies of QCD with many 



flavors [10 1 supports this result. 

The above result is quite interesting and it is desired to study the same problem by 
other methods. Recently a closely related model, the random flux model, was studied 



by a lattice formulation [21]. There a mathematical tool invented by ZirnbauerJTl] 
was used and the DoS was calculated for the single-flavor case. Then it is interesting 
to calculate the DoS for many flavors and to see if qualitative difference appears as 
the number of flavors is increased. 

In the present studies we cannot estimate the critical value of Nf for the delocalization- 
localization transition(DLT). Closely related and tractable model is a quantum U(Nf) 
ant if err magnetic Heisenberg model |E|. There a phase transition, which corresponds 
to the DLT, is predicted and the critical value of Nf is estimated. 
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